Measuring the proton spectrum in neutron decay - latest results with aSPECT 
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Abstract 

^ The retardation spectrometer aSPECT was built to measure the shape of the proton spectrum in free neutron decay with high 
precision. This aUows us to determine the antineutrino electron angular correlation coefficient a. We aim for a precision more than 
one order of magnitude better than the present best value, which is Aa/a — 5%. 



In a recent beam time performed at the Institut Laue-Langevin during April / May 2008 we reached a statistical accuracy of 
about 2 % per 24 hours measurement time. Several systematic effects were investigated experimentally. We expect the total relative 
J uncertainty to be well below 5 %. 
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• 1. Introduction 
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The decay of the free neutron allows one to determine 
the couphng constants of the weak interaction and can be 
J> p-sed to search for physics beyond the standard model of 
elementary particle physics [1,2]. The decay rate for unpo- 
larised neutrons can be described as [3] : 
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(1) 



where rric is the mass of the electron and Pc, Pu, E^ and E^, 
are the momenta and energies of electron and antineutrino, 
respectively. 

Within the standard model, where neutron decay is de- 
scribed as a V-A type interaction, the antineutrino elec- 
tron angular correlation coefficient a depends only on the 
ratio of the axial vector and vector coupling constants A = 
9a/ 9v- 
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Since A describes the renormalisation of the axial vector 
current by the structure of the nucleon, it cannot be calcu- 
lated well enough from first principles and thus has to be 
determined experimentally. 

From A and the neutron lifetime, the element Vud of the 
Cabibbo-Kobayashi-Maskawa matrix can be derived. Fur- 
thermore, the determination of A from different observables 
(e.g. from a and the beta asymmetry A [3]) permits cross- 
checks of the theory and searches for non-V-A couplings [1] . 

2. The spectrometer 

The neutrino is hard to detect, hence we infer a from the 
shape of the proton recoil spectrum. aSPECT is a retarda- 
tion spectrometer. This means, the spectrum is measured 
by counting all decay protons that overcome a potential 
barrier. By varying the height of the barrier the shape of 
the proton spectrum can be reconstructed. 

A beam of cold, unpolarised neutrons passes through the 
spectrometer, where about 10^* of the neutrons decay in 
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the decay volume (see fig.[T]). This region is held at ground 
potential, from here the decay protons are guided to the 
proton detector by a strong magnetic field (about 2T in 
the decay volume during the beam time discussed below). 
About half of the protons are emitted in the opposite direc- 
tion. These are reflected by an electrostatic mirror and thus 
finally all protons are directed towards the detector. The 
endpoint energy of the proton spectrum is about 751 eV, 
therefore a voltage of Um = 820 V at the electrostatic mir- 
ror reflects all protons. Before the protons can reach the 
detector, they have to overcome a potential barrier, which 
is maximal at the so-called analysing plane. The barrier po- 
tential is generated by a 54 cm long cylindrical electrode, 
held at a voltage Ua- In the analysing plane the magnetic 
field is about 5 times lower than in the decay volume. Pro- 
tons travel from the decay volume to the analysing plane, 
gyrating about a magnetic field line. As a result of the adi- 
abatic invariance of the magnetic moment, part of the pro- 
ton's momentum transverse to the field line is transferred 
into parallel momentum. In the adiabatic approximation, 
the probability that a proton overcomes the potential bar- 
rier can be analytically calculated as a piecewise function: 
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where T is the kinetic energy and e the charge of the proton. 
This transmission function i^tr(r) depends only on Ua and 
the ratio of the magnetic fields in the analysing plane and 
the decay volume, = Ba/Bq. 

Further details of the spectrometer design and systematic 
effects may be found in [4, 5]. 

After the electrostatic barrier the magnetic field increases 
to focus particles that overcome the barrier onto a silicon 
drift detector (SDD) [6]. The detector is held at a high neg- 
ative potential to post-accelerate the protons to detectable 
energies, and to ensure that the protons overcome the mag- 
netic mirror created by the increasing field. A SDD is a 
semiconductor detector based on the principle of sidewards 
depletion [7] , which allows the depletion of a large detec- 
tor volume with a small readout node. In the spectrometer, 
one chip with 3 detector pads (10 x 10 mm^ each) is used. 
The first amplifying transistor as well as a diode for tem- 
perature measurement are implemented on the chip PI The 
detector signal is amplified and digitised by a continuously 
sampling 12 bit analogue to digital converter with a 20 MHz 
clock. The data is then analysed on-line with respect to 
a triggering algorithm by a FPGA (Field Programmable 
Gate Array) which sends a 5 /iS long event to a storage PC 
in case of a trigger. This event length is the dead time of 
the detection system. 
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Results with a smaller test detector have been published in [8]. 



Figure 1. Simplified sketch of the aSPECT spectrometer (a detailed 
picture can be found in [5]). See text for description. 

The SDD allowed us to significantly reduce the acceler- 
ation voltage, compared to the previously used silicon PIN 
diode detector. Thus we avoided problems like electrical 
breakdowns or instabilities of the background due to field 
emission. During the beam time post-acceleration voltages 
from —10 to — 15kV were used. Some pulse-height spectra 
with different acceleration voltages are shown in fig.IH 

Protons hit the detector with slightly different incident 
angles and energies due to the distribution of initial angles 
and energies, thus an angle or energy dependence in the 
detection probability leads to a systematic effect. Although 
first simulations show that this effect is rather small, it will 
be investigated experimentally at a proton source dedicated 
for detector tests [9] . 

aSPECT was set up at the cold neutron beam line PFIB 
[10] at the Institut Laue-Langevin. The neutron beam was 
shaped by apertures made from sintered boron carbide 
in the beam tube before the spectrometer (pressure ~ 
lO^^mbar) and apertures from isotopically enriched ^LiF 
inside the spectrometer (~ 5 x 10~^ mbar). The vacua were 
separated by 250 fim thick MgA13Znl windows. A fast neu- 
tron shutter in front of the beam tube permitted to open 
and close the neutron beam during data taking. Behind the 
spectrometer, a beam stop from boron carbide with an in- 
tegrated neutron flux monitor was installed. The thermal 
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Figure 2. Typical pulse-height spectra from the beam time, taken 
with the silicon drift detector at different acceleration voltages. By 
increasing the acceleration voltage the proton peak (right peak) is 
shifted to higher ADC channels, whereas the position of the elec- 
tronics noise (left peak) is not influenced. 

neutron flux at the exit window of the spectrometer was 
6 X lO^n/cm^s. Outside the spectrometer the beam was 
shielded by several layers of ^"B loaded rubber and 10 cm 
of lead. 

3. Investigations of systematic effects 

All information in this section is preliminary. 

The proton spectrum was determined by setting the 
analysing plane electrode to 7 different voltages Ua- 50 V 
is used for the normalisation of the count rates, 400 V ap- 
proximately provides the best statistical sensitivity towards 
a [11]. The background was quantified by measurements 
with 780 V. Most of the background is caused by beam- 
related electrons from neutron decay which are also guided 
towards the detector. At 780 V no proton can overcome the 
barrier, whereas the background is only marginally influ- 
enced. The pure proton count rate is obtained by subtract- 
ing the background from the measurements with lower volt- 
ages. To learn about systematic effects and to gain a more 
precise knowledge of the shape of the spectrum V, 50 V, 
250 V, 400 V, 500 V, 600 V, and 780 V were measured. 

Typically, about 470 events per second with pulse- heights 
in the integration window from 160 to 900 ADC channels 
(fig-Hi) at 50 V analysing plane voltage were counted on 
one detector pad. At Ua = 780 V the count rate in the 
same window was about 7 counts per second. This means, 
approximately 460 protons per second were detected. With 
closed neutron shutter the count rate dropped to 0.2 counts 
per second in the integration window. 

In our previous beam time, the background count rate 
without neutron beam at Ua — 780 V was higher than at 
Ua — 50 V by several Hz [5]. For this beam time, changes 
were made to the detector high voltage and the rest of 
the electrodes system to suppress the dependence of the 
background count rate on the analysing plane voltage. So 



far, such a dependence was not clearly identified in the 
present beam time. For example, for a run of 16 hours, the 
background count rate measured during the periods with 
closed neutron shutter was 0.197(18) Hz for Ua — 50 V and 
0.154(20) Hz for Ua = 780 V. 

The absolute height of the potential barrier has to be 
known precisely to determine the transmission function, 
eq. ^ . The voltage applied to the electrode is monitored 
by a precise multimetei[fj. The accuracy of the voltage set- 
tings was better than 5 mV, limited by the multimeter's cal- 
ibration. Stability and reproducibility of the voltage were 
much better. However, the electrostatic potential inside the 
electrode is affected by the work function of the electrode 
surface. In a cylindrical sample electrode at room tempera- 
ture, a variation of up to 100 mV was found. To investigate 
this effect, further measurements with a Kelvin probe are 
ongoing. 

Protons and positive ions with too low energy to over- 
come the potential barrier are trapped between the electro- 
static mirror and the analysing plane. They are removed 
by an E X B electrode between the decay volume and the 
analysing plane (see fig. [J). This electrode consists of two 
half cyHnders set to different voltages Uexb,i and J7exb,2- 
Trapped particles oscillate between the mirror and the 
analysing plane. Every time they pass the electrical field 
produced by the E x B electrode they drift into the same 
direction, perpendicular to both E and B. 

In addition, electrons and negative ions may be trapped 
in the region of the analysing plane electrode. 

These trapping effects were investigated with two meth- 
ods: 

(i) For each analysing plane voltage Ua, the measure- 
ment started with closed neutron shutter during a 
time ti , then the shutter was opened for t2 , and closed 
again for ^3. This mode allowed us to investigate 
the environmental, not beam-related background and 
possible differences in the background count rate be- 
fore and after the shutter was opened. A higher count 
rate in the last part would point to either trapped par- 
ticles inside the spectrometer or activation of some 
material. Typically, the total measurement time with 
closed shutter was about half as long as the time with 
open shutter (ti + ^ ^'-a)- 

(ii) The analysing plane was set to 780 V so that no pro- 
tons should be able to overcome the barrier and the 
E X B drift potential was reduced from the standard 
setting of Uexbs I Uexb.2 = -1000 V | - 50 V on the 
two half cylinders of the electrode in several steps 
down to —2.5 V j V. A comparison of two measure- 
ments with -2.5 V I V and -200 V | V is shown in 
fig-El For the lowest drift voltage of 2.5 V there are 
only some counts before the shutter is opened, this is 
the environmental background. As soon as the shut- 
ter is opened, protons are trapped and the count rate 
starts to fluctuate. Even after the shutter is closed 
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Figure 3. The background count rate in the integration window for 
two different settings of the E X B electrode. The dashed hues show 
where the neutron shutter was opened (after ~ lis) and later closed 
again (~ 52 s). 

again, there are still particles that reach the detec- 
tor. Already a potential of — 10V| OV reduces this 
behaviour strongly, and with —200 V | V the count 
rate is stable as long as the shutter is open and drops 
back to the environmental background count rate im- 
mediately after the shutter is closed again. 
The knowledge of the magnetic field in both decay vol- 
ume and analysing plane is important for the calculation of 
the transmission function, cq. ([3|). To determine the mag- 
netic field ratio re, the magnetic field has been measured 
in the spectrometer installed at the beam position, before 
and after the beam time. There was no significant change 
in the ratio or the absolute values of the fields. To check 
the response of the spectrometer, the ratio was changed by 
about 1 % by additional, external coils. These coils do not 
change the homogeneity or the shape of the magnetic field 
in both analysing plane and decay volume. As expected, the 
count rates for the different measured voltages Ua changed 
due to the altered transmission function, but the value for 
a is, within statistical errors, still the same if we use the 
adapted ratio of the magnetic fields Tb in eq. ([3]). 

To monitor the ratio re , it is planned to install a nuclear 
magnetic resonance system with polarised "^He. The NMR 
frequency will be measured in two cells close to decay vol- 
ume and analysing plane, respectively. The ratio of those 
frequencies is proportional to tb . 

The superconducting magnet coils show some hysteresis 
which might alter the field values if the magnet is set to a 
different strength. Measurements showed that despite the 
hysteresis magnetic field settings can be reproduced if for 
each change in the coil setting the same procedure is used, 
which starts with warming up the coils above their critical 
temperature. Then they are cooled down again and finally 
ramped directly to the desired current. 

Another systematic effect we investigated is the so-called 
edge effect. The neutron beam profile in the decay volume 



is projected onto the detector by the magnetic field. Close 
to the edge of the detector, two cases may occur: On the 
one hand, protons produced outside the directly projected 
volume can hit the detector, since they gyrate about a 
magnetic field line. This increases the detected count rate. 
On the other hand, protons produced within the projected 
area might miss the detector, resulting in a loss of detected 
count rate. As the gyration radius of a proton depends on 
its momentum, the proton count rate gain and loss is mo- 
mentum dependent. As long as the beam profile is perfectly 
flat and the detection efficiency of the proton detector is 
uniform the two effects cancel. But if the profile is non-fiat 
the probabilities for the two cases become different. This 
leads to an energy-dependent systematic effect. 

The effect can be calculated using Monte Carlo simula- 
tions. Still, the neutron apertures were simulated and de- 
signed to obtain a beam profile as flat as possible to min- 
imise the effect. To verify our calculations of the edge ef- 
fect, data was taken for different beam profiles. These pro- 
files were generated with an additional aperture in front of 
the spectrometer. Each profile was measured by copper foil 
activations in front of and behind the spectrometer. 

4. Conclusion 

In the recent beam time we showed that our spectrome- 
ter is capable of measuring the proton recoil spectrum with 
high precision. Compared to a previous beam time, the 
major improvements were the new proton detector, the re- 
design of several electrodes, and a better vacuum. From the 
investigated systematic effects and the collected statistics, 
we expect a total relative error well below 5 %. 

With the knowledge of the systematic effects gained in 
this beam time we should be able to improve aSPECT to 
permit a measurement of a with its design accuracy, which 
is 0.3 % in a further beam time. 
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